Acceleration and vortex filaments in turbulence 
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We report recent results from a high resolution numerical study of fluid particles 
transported by a fully developed turbulent flow. Single particle trajectories were 
followed for a time range spanning more than three decades, from less than a tenth 
of the Kolmogorov time-scale up to one large-eddy turnover time. We present some 
results concerning acceleration statistics and the statistics of trapping by vortex 
filaments. 
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Lagrangian statistics of particles advected by a turbulent velocity field, u{x, t), are impor- 
tant both for their theoretical implications (l| and for applications, such as the development 
of phenomenological and stochastic models for turbulent mixing yl. Despite recent advances 
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in experimental techniques for measuring Lagrangian turbulent statistics p|, y, y, la| , direct 
numerical simulations (DNS) still offer higher accuracy albeit at a slightly lower Reynolds 
number BBS. Here, we describe Lagrangian statistics of velocity and acceleration in 
terms of the multifractal formalism. At variance with other descriptions based on equi- 
librium statistics (see e.g. Q, H Q, critically reviewed in jl3|), this approach has the 
advantage of being founded on solid phenomenological grounds. Hence, we propose a deriva- 
tion of the Lagrangian statistics directly from the Eulerian statistics. 

We analyze Lagrangian data obtained from a recent Direct Numerical Simulation (DNS) 
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of forced homogeneous isotropic turbulence |l4l.llq which was performed on 512^ and 1024^ 
cubic lattices with Reynolds numbers up to Rx ~ 280. The Navier-Stokes equations were 
integrated using fully de- aliased pseudo-spectral methods for a total time T ^ Tl- Two mil- 
lions of Lagrangian particles (passive tracers) were injected into the flow once a statistically 
stationary velocity field had been obtained. The positions and velocities of the particles were 
stored at a sampling rate of 0.07r^. The velocity of the Lagrangian particles was calculated 
using linear interpolation. Acceleration was calculated both as the derivative of the particle 
velocity and by direct computation from all three forces acting on the particle (i.e. pressure 
gradients, viscous forces and large scale forcing): the two measurements were found to be 
in very good agreement. Finally, the flow was forced by keeping the total energy constant 
in each of the first two wavenumber shells. For more details on the simulation, see 3| ■ 



I. VELOCITY AND ACCELERATION STATISTICS 

Velocity statistics along a particle trajectory can be measured by means of the Lagrangian 
structure functions, >S'p(r) = (((5^t>)P) where 6rV is the Lagrangian increment of one compo- 
nent of the velocity field in a time lag r. A simple way to link the Lagrangian velocity 
increment, 6rV, to the Eulerian one, 6rU, is to consider the velocity fluctuations along a 
particle trajectory as the superposition of different contributions from eddies of all sizes. In 
a time-lag r the contributions from eddies smaller than a given scale, r, are uncorrelated, 
and we may write 6rV ~ 6rU. Assuming that typical eddy turn over time r at a given spatial 
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scale r can be expressed as r^. ~ r/SrU, one obtains: 

rh 

6rV r^drU r ~ ^r^-^ (1) 

where h is the local scaling exponent characterizing the Eulerian fluctuation in the mul- 
tifractal phenomenology 16|. Also, Lq^vq are the integral scale and the typical velocity, 
respectively. With respect to the the usual multifractal phenomenology of fully developed 
turbulence, the presence of a fluctuating eddy turn over time is the only extra additional 
ingredient to take into account in the Lagrangian reference frame. 

Using (0), one can estimate the Lagrangian velocity structure function: 

hp+Z-D{h) 

WM<)/^,dA(^) *-* . (2) 

where the factor (r/T/^)^^^-^^^))/^^^^) is the probability of observing an exponent /z. in a time- 
lag r, and D{h) is the dimension of the fractal set where the exponent h is observed. The 
Lagrangian scaling exponents Chip) can be estimated by a saddle point approximation, for 

a(.)^.f(^-i^). (3) 

We would like to stress that for the D{h) curve we have chosen that of the Eulerian statis- 
tics. In other words, the prediction © is free of any additional parameter once the Eulerian 
statistics are assumed Il4 llTI. llSll. 
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In Fig. (P), we present the Extended Self Similarity (ESS) [19| log-log plot of 5'p(r) versus 
5*2 (r) as calculated from our DNS. The logarithmic local slopes shown in the inset display 
a deterioration of scaling quality for small times. We explain this strong bottleneck for 
time lags, r G [r^, lOr^], in terms of trapping events inside vortical structures [3]: a dy- 
namical effect which may strongly affect scaling properties and which a simple multifractal 
model cannot capture. For this reason, scaling properties are recovered only using ESS 
and for large time lags, r > lOr^. In this interval a satisfactory agreement with the mul- 
tifractal prediction Q is observed, namely from the multifractal model one can estimate 
Cl(4)/Cl(2) = 1.71, Cl(6)/Cl(2) = 2.26, Cl(8)/Cl(2) = 2.72 while from our DNS we measured 
a(4)/CL(2) = 1.7 ± 0.05, Cl(6)/Cl(2) = 2.2 ± 0.07, a(8)/CL(2) = 2.75 ± 0.1. 

A similar phenomenological argument can be used to make a prediction for the accelera- 
tion probability density function (pdf). The acceleration can be deflned as: 



(4) 
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FIG. 1: ESS plot of Lagrangian velocity structure function S'p(t) versus S2{t). Symbols refer 
to the DNS data for p = 8,6,4 from top to bottom. Lines have slopes Cl(p)/Cl(2) given by 
the multifractal prediction ^ with a D{h) curve taken from the She-Leveque prediction ^20]. In 
the inset, we show the local slopes versus time t/t^, and their comparison with the respective 
multifractal prediction (straight lines). 

As the Kolmogorov scale itself, rj, fluctuates in the multifractal formalism: ri{h,V()) ~ 
(uLq/voJ \ so does the Kolmogorov time scale, Trj{h,vo). Using (jT)) and (jH) evalu- 
ated at ?7, we get for a given h and Vq'- 
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a{h, vo) ~ u—v^^'' Lo . (5) 

The pdf of the acceleration can be derived by integrating (jSJ over all h and vq, weighted 
with their respective probabilities, {T^{h,vo)/TL{vQ))^'^^^^^^^^^^~^^ and V{vo). It remains to 
specify a form for the large scale velocity pdf, which we assume to be Gaussian: V{vq) = 
1/ \J 271(7'^ exp{—VQ / 2crl) , where = (vq). Integration over vq gives: 

Via) / aha s u 3 Lq^ ' cr,, x 

Jh£l 

In order to compare the DNS data with the multifractal prediction we normalize the acceler- 
ation by the rms acceleration a a = {a^Y^"^ oc R\. In terms of the dimensionless acceleration. 
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FIG. 2: Log-linear plot of the acceleration pdf. The crosses are the DNS data, the solid black line 
is the multifractal prediction and the green line is the K41 prediction. The statistical uncertainty 
in the pdf is quantified by assuming that fluctuations grow proportional to the square root of the 
number of events. Inset: a^V{a) for the DNS data (crosses) and the multifractal prediction. 
a = ajoai (0) becomes 

/■ , (fe-5+-D(ft)) (fe) / 1,2(1+;.) (feA 

V[a) a 3 ^exp — -a a R^^ '] ah, (7) 

where = xih-5 + D{h)) /6 + 2{2D{h) +2h-7) /3, z{h) = x(l + /i)/3 + 4(2/i - l)/3 and 
X = sup^j {2{D{h) — Ah — + h) ) . For more details on how the numerical integration of 
© is made we refer the reader to jlq]. 

In Fig. (j2|) we compare the acceleration pdf computed from the DNS data with the 
multifractal prediction (jZj). The large number of Lagrangian particles used in the DNS 
(~ 10^) allows us to detect events up to SOctq. The accuracy of the statistics is improved by 
averaging over the total duration of the simulation and all spatial directions, since the flow 
is stationary and isotropic at small-scales. Also shown in Fig. Q is the K41 prediction for 
the acceleration pdf VxAiio) ~ a~^/^exp (^—a^/^/2}j which can be recovered from ((7j) with 
h = 1/3, D{h) = 3 and x = 1- As evident from Fig. (j2)), the multifractal prediction ((Zj) 
captures the shape of the acceleration pdf much better than the K41 prediction. What is 
remarkable is that (0) agrees with the DNS data well into the tails of the distribution - from 
the order of one standard deviation cTq up to order 70aa- This result is obtained using the 
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FIG. 3: (Left panel) Trajectories with an intense value of the acceleration have been selected: as 
it can be seen, this corresponds to select tracers trapped into vortex filaments. Arrows and colors 
encode the velocity (magnitude and direction) of the particle. Rendering is realized with OpenDX. 
The movie (see multimedia enhancement) shows the flow as seen by riding this particle, before and 
during the trapping event. (Right panel) We show, in natural units, the behavior of one component 
of the centripetal and of the longitudinal acceleration (for details see the text). Notice the strong 
sign persistence of the centripetal acceleration with respect to the longitudinal one. 

She-Leveque model for the curve D{h) 




II. ACCELERATION TAILS AND SPIRALING MOTION 

This and previous work j^, 0, Q| has collected evidence which highlights the relevance 
to Lagrangian turbulence of strong spiraling motions corresponding to trapping events, i.e. 
passive particles trapped in small scale vortex filaments. So we identify the strong bottleneck 
effect visible in Figure Q and also the presence of extremely rare fluctuations in the pdf of 
the acceleration (see Figure |21). To illustrate better these strong events, we plot one of them 
in Figure 121 As is evident, the particle while moving slowly and smoothly, at some point gets 
trapped in a vortex filament and starts a spiraling motion characterized by huge values of 
the acceleration and by a "quasi-monochromatic" signal on all the velocity field components. 
Here, we suggest a way to characterize such events. This is of course a difficult task because 
not all the "trapping events" are so clearly detectable as that shown in Figure El 
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Indeed the motion of a particle in a turbulent field will be characterized by different 
accelerations and decelerations, not necessarily associated with spiraling motion (on average 
the mean value of the acceleration will be zero). In a spiraling motion the velocity v and 
acceleration a are orthogonal. Furthermore in a circular uniform motion the angular velocity, 
u, can be related to the centripetal acceleration Oc = w^r and to the linear velocity v = ur. 
We expect that in trapping events such as the one depicted in Fig. Q the centripetal 
acceleration is intense and much more persistent than the longitudinal acceleration (i.e. the 
acceleration in the direction of the motion). To make this statement quantitative, we have 
studied the average of the centripetal, ac = axv = ax-^^, and longitudinal acceleration, 
ai = {a ■ v)v, over a time window which can vary up to 9r^, A = {0.1, 3, 9}r^: 

af(t) = (a,)^ = i^*^^dt'a,(t'); (8) 
af (t) = (a,)^ = - dt'a,{t'). (9) 

We expect that the pdfs of the averaged centripetal and longitudinal acceleration will behave 
very differently with increasing the window size, A. In particular, the strong persistence 
of the centripetal acceleration up to lOr^ suggests that the centripetal pdf V{a^) should 
remain almost unchanged when varying A, while the longitudinal one V{af^) should become 
less and less intermittent. This is what we show in Fig. 

In order to investigate further the role of trapping in vortices, we can define a typical 
radius of gyration Tc and its typical eddy turnover time Tc, as: 

112 II 

\v\ \v\ 

fc = 1 TT and Tc = — (10) 

\a X v\ \a X v\ 

Notice that using axv corresponds to selecting the centripetal values of the acceleration and 

hence augmenting the signal/noise ratio of spiraling motions with respect to the background 

of turbulent motions. The previous expressions applied to a typical vortex filament give 

Tc ~ and ~ r^. Similarly one may define a typical time based on the "longitudinal 

acceleration" : ti = |i;|/|(a ■ v)v\. Incoherent fluctuations with typical times of the order of 

should be averaged out once we measure the mean centripetal and longitudinal accelerations 

averaged over a window with A > r in expression (jlO|) . On the other hand, the signal coming 

from coherent vortex should not be affected by the averaging procedure and keeps its value: 

as a consequence, we should see events with Tc ~ even upon averaging. Going through 

Figure El we can observe, with increasing window size, the different behaviors of the pdfs of 
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FIG. 4: Pdf of the averaged centripetal Oc (left panel), and longitudinal ai (right panel) acceleration 
components. The acceleration is averaged over a time window of size A = {0.1, 3, 9}r^ (respectively 
corresponding to colors red, green and blue). 




FIG. 5: Pdf of the characteristic time estimated from the centripetal (in red) and longitudinal 
(in green) accelerations (in units of r^) V{tc) and V{ti) respectively, for A = O.lr^ (left panel); 
A = 3r^ (central panel); A = 9r^ (right panel). 

the centripetal and longitudinal characteristic times, Tc and r/ respectively. It is interesting 
to notice that the left tail of the centripetal pdf is quite robust, showing the presence of 
characteristic times of the order of ~ even after averaging over a window with A = 9r^. 
On the other hand the longitudinal characteristic times of order ti ~ soon disappear as 
long as A > r^. We interpret this as further evidence of the importance of trapping in vortex 
filaments. 
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III. CONCLUSIONS 

We have presented results on the Lagrangian single-particle statistics from DNS of fully 
developed turbulence. In particular we have shown that (i) in the large time lag limit, 
IOtjj < t < Tl, velocity structure functions are well reproduced by a standard adaptation 
of the Eulerian multifractal formalism to the Lagrangian framework; (ii) the acceleration 
statistics are also well captured by the multifractal prediction; (iii) for time lags of the order 
of the Kolmogorov time scale, r^, up to time lags lOr^, the trapping by persistent vortex 
filaments may strongly affect the particle statistics. The last statement is supported both 
by the scaling of the Lagrangian statistics and by a new analysis based on the centripetal 
and longitudinal acceleration statistics. 
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